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We study the contribution to the decay of de Sitter vacua from bounces with 
0(3) x 0(2) symmetry. These correspond to the thermal production of a vacuum 
bubble at the center of a horizon volume with radius rjj and a temperature defined 
by the horizon. They are analogues of the flat spacetime bounces, independent 
of Euclidean time, that correspond to thermal production of a critical bubble. If 
either the strength of gravity or the false vacuum energy are increased, with all 
other parameters held fixed, the bounces approach, and eventually merge with, the 
Hawking-Moss solution. Increasing the height of the barrier separating the true and 
false vacuum, and thus the tension in the bubble wall, causes the center of the bubble 
wall to approach, but never reach, the horizon. This is in contrast with the prediction 
of the thin-wall approximation, which inevitably breaks down when the wall is near 
the horizon. Our numerical results show that the Euclidean action of our solutions is 
always greater than that of the corresponding 0(4)-symmetric Coleman-De Luccia 
bounce. 
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I. INTRODUCTION 

Vacuum decay by quantum tunneling in curved spacetime was first studied more than 
thirty years ago by Coleman and De Luccia (CDL) [l], using a formalism based on bounce 
solutions of the Euclidean field equations. In recent years there has been renewed interest 
in this topic, motivated in part by its relevance for a possible string landscape. This has 
uncovered new aspects of the problem that become significant if the mass scales approach 
the Planck mass or the bounce radius approaches the horizon length. 

It was assumed in pj that the dominant bounce solution had an 0(4) symmetry; we will 
refer to such solutions as CDL bounces. In this paper we will investigate a different class of 
solutions, with an 0(3) x 0(2) symmetry, that can occur when the tunneling is from a de 
Sitter spacetime. 

Our motivation for studying these is two-fold. For the case of tunneling in flat spacetime 
at zero temperature it has been shown that the bounce with minimum action in a theory 
with a single scalar field has an 0(4) symmetry, and it is generally believed that this result 
extends to theories involving more fields. Indeed, it is quite possible that all bounce solutions 
in these theories are 0(4)-symmetric. No proof of a comparable result in curved spacetime 
has ever been obtained. The 0(3) x 0(2) symmetry that we impose (which would not be 
possible in flat spacetime at zero temperature) reduces the field equations for the bounce 
to a set of ordinary differential equations, thus making the problem quite tractable and 
allowing comparison of the actions of the two types of bounces. 

A second motivation is that these solutions shed light on a reinterpretation of the CDL 
prescription that was advocated recently by Brown and one of us (2|. It was argued there 
that in the case of decay from a de Sitter vacuum the bounce should be understood as 
describing tunneling in a finite horizon volume at finite temperature. This has implications 
for extracting from the bounce the configuration immediately after tunneling. The CDL 
bounce is topologically a four-sphere, and it is commonly asserted that the configuration 
after tunneling is obtained from the three-sphere lying on the "equator" of the bounce. 
By contrast, it was argued in Ref. [2| that half of this three- sphere gives the configuration 
within the horizon volume after tunneling, with the other half giving the configuration before 
tunneling. The fields beyond the horizon, either before or after tunneling, are never specified 
and do not enter the calculation. For CDL bounces corresponding to vacuum bubbles much 
less than horizon size there is little practical difference between the two interpretations. The 
differences become significant as the bubble approaches horizon size and, especially, in the 
Hawking-Moss (if limit. For the bounces we consider, the differences in interpretation are 
manifest for all bubble sizes. From the viewpoint of Ref. j2[ these bounces correspond to 
thermal production of a single bubble, whereas the conventional viewpoint would yield the 
somewhat surprising conclusion that a pair of bubbles had been produced. 

Our interest in this these bounces was first raised by the work of Garriga and Megevand [H, 
. Motivated by the analogy with the thermal production of bubbles in flat spacetime, they 
studied bounces with 0(3) x 0(2) symmetry and, indeed, did interpret these as corresponding 
to the nucleation of a pair of bubbles. The discussion was originally in the context of brane 
nucleation 0, 0| , in which the bubble walls are closed 2-branes, but their analysis can also be 
applied to field theory nucleation in the thin-wall approximation. We will differ from them 
in not restricting ourselves to this thin- wall limit. In fact, we will see that the thin- wall 
approximation inevitably breaks down when the bubble becomes sufficiently large, in what 
is perhaps the most interesting limit of these solutions. 
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We work in the context of a scalar field theory with Lagrangian 

£=^0) 2 -V(0), (l.l) 

where V(<ft) has two local minima, with the higher (lower) being termed the false (true) 
vacuum. We assume that both minima have positive energy, and are thus de Sitter space- 
times. Although our focus is on the nucleation of a true vacuum bubble in a false vacuum 
background, we also consider the nucleation of a bubble of false vacuum within a de Sitter 
true vacuum background. Although this process is energetically forbidden in the absence of 
gravity, it can occur in de Sitter spacetime with its nonzero temperature T^s |8(. In contrast 
with the CDL case, the bounces governing the nucleation of true vacuum and false vacuum 
bubbles are not identical. 

The remainder of this paper is organized as follows. In Sec. II, we review the essentials of 
vacuum tunneling, focusing on some features of particular relevance here. In Sec. Ill, we set 
up the formalism that we will use, and obtain some basic results concerning our solutions. In 
Sec. IV we discuss the application of the thin-wall approximation to our problem, explaining 
how its applicability is more limited than in the case of bounces with 0(4) symmetry. We 
then turn to the description of the 0(3) x 0(2)-symmetric bounces themselves. Some limiting 
cases that can be handled analytically are discussed in Sec. V, and in Sec. VI we describe 
the results of our numerical solutions. Section VII contains some concluding remarks. 



II. REVIEW OF VACUUM TUNNELING 

A. Flat spacetime at zero temperature 

In flat spacetime at zero temperature, vacuum decay occurs through the production 
of bubbles of true vacuum via quantum tunneling. Importantly, this is not a tunneling 
process from homogeneous false vacuum to homogeneous true vacuum, but rather tunneling 
from an initial configuration that is pure false vacuum to a final configuration containing 
a bubble of true vacuum (or an approximation thereof) within a false vacuum background. 
The tunneling is not actually through the one-dimensional barrier that appears in a plot of 
V ((/>), but rather through an infinite-dimensional barrier in configuration space defined by 
the potential energy functional 



u[4>(x)] = J 



d 3 x 



i(V M 0) 2 + V(0) 



(2.1) 



One approach to this problem is via the WKB method. With one degree of freedom, 
q(t), the WKB approximation gives an amplitude for tunneling through a barrier that is 
proportional to e~ B , where B is of the form 

B=\r dqp[V{q)-E] (2.2) 

with V(qi) = V(q2) = E. For a barrier in a multidimensional configuration space, one 
must find a path through the barrier for which the corresponding one-dimensional tunneling 
exponent B is a local minimum j§, 10]. Coleman [ll[ showed that the problem of finding a 



stationary point of B is equivalent to finding a "bounce" solution of the Euclidean equations 
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FIG. 1: Schematic view of a bounce in flat spacetime. The shaded area is the region of approxi- 
mately true vacuum, while the heavy circle surrounding it represents the wall region. In actuality, 
the field is not strictly constant outside the wall, but instead has an exponential tail and only 
reaches its false vacuum value at infinite distance. The horizontal dashed lines represent a se- 
quence of spatial hypersurfaces of constant r that trace out a path in configuration space. The 
lower heavy dashed line (corresponding to r = — oo) indicates the initial pure false vacuum state, 
while the heavy dashed line through the center of the bounce corresponds to the field configuration 
at bubble nucleation. 



of motion. In the field theory case, the bounce is a solution 0(x, r) in four- dimensional Eu- 
clidean space. The hypersurface at the initial value Tinit (which is equal to — oo for tunneling 
from the false vacuum) gives the configuration before tunneling, while the configuration 
0(x, r = 0) along the hypersurface through the center of the bounce gives the optimal exit 
point from the potential energy barrier; i.e., the form of the bubble at nucleation. The 
hypersurfaces at intermediate Euclidean times — oo < r < give successive configurations 
along the optimal tunneling path. Those for < t < oo are simply mirror images of those 
for negative r, thus motivating the term "bounce". The bubble nucleation rate per unit 
volume is of the form F = Ae~ B , with B being equal to the difference of the Euclidean 
actions of the bounce and the pure false vacuum. 

Subject to only mild assumptions on the potential [l2|, the bounce of minimum action 
can be shown to have an 0(4) symmetry. This bounce has a central region, in which <ft is 
close to its true vacuum value, within an infinite region of false vacuum, with the two regions 
separated by a wall of finite thickness in which traverses the potential barrier, as shown 
in Fig. [TJ 

An alternative approach is to use a Euclidean path integral to calculate the imaginary 



part of the energy of the false vacuum, which is closely related to T [13j. This gives the 
prefactor A in terms of the eigenvalues of small fluctuations about the bounce, with the 
imaginary part of the energy emerging from the presence of a single negative eigenvalue. If 
there are additional negative eigenvalues, the path specified by the bounce is not even a local 
minimum of £>, but rather a saddle point, and the bounce solution must be discarded [14]. 
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B. Finite temperature 



Both the WKB and the path integral approaches show that tunneling at a nonzero tem- 
perature T is governed by a bounce solution that is periodic in imaginary time with period 
P = 1/T. Two types of periodic bounces can be distinguished. 

At low temperature the bounce is a somewhat deformed and elongated version of the 
zero-temperature bounce, with the 0(4) symmetry reduced to an 0(3) symmetry. Such 
bounces correspond to thermally assisted tunneling processes in which thermal effects allow 
the system to start the tunneling at a somewhat higher point on the potential energy barrier. 
At sufficiently high temperatures this solution may cease to exist. 

The second type of solution corresponds to a thermal fluctuation to a saddle point at 
the top of a "mountain pass" over the potential energy barrier. These solutions are in- 
dependent of r. Each three-dimensional spatial slice is infinite in extent and describes an 
0(3)-symmetric configuration 0(x) with a central true vacuum bubble of critical size. The 
Euclidean action of these solutions is equal to the energy of this critical bubble divided by 
the temperature. Although these solutions exist for all values of T, at low temperature they 
have multiple negative modes and must therefore be discarded. Because this point will be 
important for us later, let us examine it a bit more closely. 

A critical bubble solution 0(x) is a saddle point of the potential energy functional U[(f>]. 
The second variation of U has a single negative mode ?/>-(x) satisfying 

V- = -V 2 V-(x) + \/"(0(x))^(x) = A_V-« (2.3) 
with A_ = — k 2 < 0. The second variation of the Euclidean action, 

f = -£-V + V W x)), (2.4) 
must then have a spectrum that includes eigenmodes of the form 

cos(^)iMx) n = 0,l,... (2.5) 

and 

sin(^^-(x) n = l,2,... (2.6) 



with eigenvalues 



A„ = W^Y. (2.7) 



Hence, at temperatures T < k/2n there are multiple negative modes, and the r-independent 
bounce must be discarded. In the flat spacetime thin-wall approximation with a critical 
bubble of radius R one finds that k = y/2/R. More generally, we expect that k ~ 1/TZ, 
where 1Z is a characteristic length scale of the critical bubble solution. Thus, we must require 
that 

^sk- < 2 - 8 > 
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FIG. 2: Two views of the CDL bounce. The bounce has the topology of a four-sphere, and is shown 
here as if embedded in a five-dimensional Euclidean space. Because of the bounce's 0(4) symmetry, 
the fields and metric have nontrivial dependence only on x^. The left-hand view shows the bounce 
as viewed from a point outside the bounce on an axis perpendicular to the X4-X5 plane. The 
shaded area corresponds to the region of approximately true vacuum. The horizontal dashed line 
has often been taken, by analogy with the heavy dashed line in FigfTJ as giving the configuration 
at bubble nucleation. The right-hand view shows the bounce as viewed from a point to its right on 
the positive xs-axis. A view from a corresponding point on the negative 25-axis would show only 
a region of approximate false vacuum. 



C. Including gravity 



Coleman and De Luccia proposed jjLf| that the effects of gravity on bubble nucleation could 
be obtained by including a curvature term in the Euclidean action and seeking a bounce 
solution that satisfies the coupled Euclidean Einstein and matter equations. If one assumes 
that the bounce has an 0(4) symmetry, and V(<p) is everywhere positive, the bounce solution 
is topologically a four-sphere that can be viewed as the hypersurface 

x\ + x\ + x\ + x\ + x\ = H~ 2 (2.9) 

in a five-dimensional Euclidean space, although not necessarily with the metric implied 
by the embedding. An 0(4)-symmetric solution is obtained by requiring invariance under 
rotations among the first four coordinates. The scalar field <fr is then a function only of x§. 

Figure [2] illustrates such a solution for the case where the bubble nucleates with a size 
much smaller than the false vacuum horizon length . Analogy with the flat spacetime 
case illustrated in Fig. [T] might suggest that the configuration after nucleation is given by 
a slice along the "equator" of this hypersphere that cuts through the center of the true 
vacuum region, as shown in the figure, with the tunneling path given by a sequence of 
parallel horizontal slices. This interpretation becomes problematic when the bubble size 
becomes comparable to the horizon length and there is no plausible candidate for the initial 
configuration. The is especially so in the limit where the bounce becomes the homogeneous 
Hawking-Moss solution [3| and a straightforward application of the interpretation would 
imply that the bounce corresponded to a transition in which all of space fluctuated to top , 
the value at the top of the barrier in V(<p). 

In Ref. 0] it was argued that one should instead interpret the Coleman-De Luccia formal- 
ism as describing thermal tunneling within a horizon volume, with the temperature being 



(a) 



(b) 



FIG. 3: The CDL bounce viewed as low-temperature thermal tunneling within a horizon volume. 
In (a) the horizontal lines represent three-balls with their center on the left and their outer edge, 
at the horizon, on the right. The horizontal line at the bottom of the figure corresponds to 
the configuration before tunneling, while the one at the center of the figure corresponds to the 
configuration after tunneling. As in Fig. [H the intermediate horizontal lines define the optimum 
tunneling path through configuration space. In (b) the same CDL solution is represented, but 
with the path through configuration space defined by a series of radial lines. Again, each of these 
represents a three-ball, with center of the ball at the edge of the figure and the two-sphere boundary 
at the horizon at the center of the figure. The initial and final configurations correspond to the 
left solid and right dashed horizontal radii, respectively. 

the one denned by the horizon. From this viewpoint, the solutions are analogous to those 
in flat spacetime at finite temperature, but with the crucial difference that the spatial slices 
of constant Euclidean time are finite in extent. Thus, these solutions exist on a space that 
is topologically the product of a three-ball (the horizon volume) and a circle (the periodic 
Euclidean time r). This suggests a "low-temperature" bounce solution such as that shown 
in Fig. |2K The boundary of this space is a family of two-spheres, each labeled by a value of 

T. 

At first sight, this seems nothing like the CDL bounce. However, this can also be rep- 
resented as in Fig. [3Jd, with the hypersurfaces of constant r appearing as diagonal sections 
cutting halfway through a four- sphere. On each section the point at the edge of the dia- 
gram corresponds to the center of the horizon volume, while the bounding horizon radius 
two-sphere appears at the center of the diagram, x 4 = x 5 = 0. Drawn in this manner, the 
bounce looks like the four-sphere bounce of Fig. It actually is topologically a four-sphere 
if the two-sphere boundaries of the sections are all identified. 

There appears to be no a priori requirement that these two-spheres be identified, nor that 
the metric should be smooth there if they are identified. Nevertheless, the 0(4)-symmetric 
CDL bounce, with an everywhere smooth metric, can be obtained in this fashion. However, 
we see that the equatorial slice is now composed of two parts. One half gives the field 
configuration at the time of bubble nucleation, but only within the horizon volume. The 
other half gives the field configuration in this horizon volume before the tunneling process. 
Note that the bounce gives us no information about the field or metric beyond the horizon, 
other than the fact that they effectively create a thermal bath for the horizon volume. 

With this understanding of the CDL bounce as describing low-temperature thermally 
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FIG. 4: Two views of the 0(3) xO(2)-symmetric bounce. These are from the same viewpoints as 
those of the CDL bounce in Fig. [2j (The view on the right would be the same if viewed from 
along any axis in the X4-X5 plane.) On the left-hand view the radial dashed lines correspond to 
the similar lines in Fig. [3b. The bounce is independent of the Euclidean time, and the three- 
dimensional configuration corresponding to any one of these lines contains a single true vacuum 
critical bubble in the center of a false vacuum horizon volume. 

assisted tunneling, we see that the gravitational analogue of the flat spacetime high- 
temperature bounce is one that is independent of r, with the slices of constant r describing 
the critical bubble configuration that is reached by thermal fluctuation. The Hawking-Moss 
solution is the simplest example, and corresponds to a spatially homogeneous critical "bub- 
ble" that fills the horizon volume. However, there can also be critical bubbles that are 
smaller than horizon size. One such is obtained simply by rotating a CDL bounce so that 
its symmetry axis is, say, the along the X\ axis. This rotated bounce corresponds to thermal 
nucleation of a critical bubble centered at the horizon (2)]. 

We will proceed by seeking a different, and perhaps more obvious, type of bounce, namely 
one that is independent of r and spherically symmetric on each constant-r slice. This would 
correspond to thermal production of a critical bubble located at the center of the horizon 
volume. Drawn as a four-sphere, the bounce would appear as in Fig. |H We see here a 
case where the analysis of Ref. |2j leads to a prediction that is clearly different from the 
conventional interpretation. Whereas the latter would view such bounces as describing the 
production of a pair of bubbles, under the interpretation of Ref. [2] only a single bubble is 
produced. 

III. FORMALISM AND SOME BASIC RESULTS 

We are thus led to seek solutions that have nontrivial behavior in three spatial dimensions, 
but that are constant in Euclidean time. If we impose spherical symmetry on the spatial 
slices, the metric on the four- dimensional space can be written in the form 

ds 2 = B(r)dr 2 + A(r)dr 2 + r 2 dtt 2 (3.1) 

and the scalar field is a function only of r. Each spatial slice is a three-ball bounded by a 
two-sphere at r = r#, with the horizon radius r# corresponding to a zero of g rr = 1/A. In 
order that the metric be nonsingular at r = 0, we must require that A(0) = 1. We do not 
fix .B(O) directly, but instead require that it be such that the product AB is equal to unity 
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at the horizon. Generically, there is a region of approximate true vacuum centered around 
the origin, r = 0, an outer region of approximate false vacuum extending out to the horizon, 
and a wall region separating the two. 

In the analogous flat spacetime case the temperature, and hence the periodicity, can be 
chosen arbitrarily, but in our curved spacetime problem we expect the temperature to be 
specified by the metric. We therefore assume that the periodicity of the Euclidean time 
is determined by the surface gravity at the horizon of the Lorentzian counterpart of our 
Euclidean metric; i.e., 



ft 



1 

T 



An 



B' 



-1/2 



in 



(3.2) 



where a prime indicates differentiation with respect to r and the subscript indicates that all 
quantities are to be evaluated at the horizon. 1 Now note that 



dr 



(AB Y (I) + (^) (I 



(3.3) 



Because A 1 = and AB = 1 at the horizon, we have B' = (A 1 )' at the horizon and 

P = -4n[B'(r H )]-\ (3.4) 

with the minus sign arising because B' is negative at the horizon. With a uniform vacuum 
ft would be 1nru\ for the bounce solutions it remains of this order of magnitude. 

With this periodicity, the four-sphere obtained by identifying all of the boundary two- 
spheres is a smooth manifold, just as in the CDL solution. However, as in the CDL case, 
we view this as incidental and do not require to be smooth on this two-sphere. Such 
smoothness would require that 0' = d<p/dr vanish at r = r#. As we will see, this is only 
possible in the (unattainable) limit in which the thin- wall approximation is exact. Although 
having </>'(?"#) ^ implies a discontinuity in slope along a path passing through the two- 
sphere, this discontinuity does not imply any cost in action, and we see no problem with 
this. 

The Euclidean action is 



S 



d 4 Xy/g 



1 



16nG 



R + C 



4vr/3 / drr 2 VAB 
Jo 



matter 
1 



lQnG 2 



k'\2 



A 



The curvature scalar that follows from the metric of Eq. ( 13.1 p is 



R 



d 



^y/ABdr 



B' 



Vab 



2 d 
r 2 dr 



(3.5) 



(3.6) 



Substituting this back into the action gives 



c_ ft 



B' 



Vab 



r=rn 



r=0 



+4nf3 f r " dr VAB 
Jo 



1 d 



8nGdr 



+ r A 



2 A 



t\2 



+ V(<f>) 



(3.7) 



Note that /3 only enters our calculations as an overall multiplicative factor in the action, and does not 
otherwise affect the solutions. It would be a straightforward matter to extend our results to an arbitrary 
temperature, although the physical origin of the temperature would be less compelling. 
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Although this is written in terms of the three functions A(r), B(r), and 0(r), it is more 
convenient to treat AB, rather than B, as the independent variable. Varying the action 
with respect to the independent variables then leads to 2 







1 d 
8^G~d^ 







"1 




+ r 2 








2 



Ll\2 



= ™.-l*GrW 







Vab 

1 d 



r 2 VAB^ 



2 A 



V(0) 



(3.8) 
(3.9) 
(3.10) 



For any solution of these equations [or even any configuration satisfying Eq. ( 13. 8p ], the 
integral in Eq. ( 13.71) vanishes. Evaluating the first term with the aid of Eq. ( 13. 41) . and recalling 
our convention that AB = 1 at the horizon, we obtain 



S ] 



bounce 



7T 

— r 
G 



H 



(3.11) 



This agrees with the thin- wall result of Ref . [5] , and could in fact have been anticipated on 
more general grounds 15|, [16 . 

Finally, the exponent in the tunneling rate is the difference between the Euclidean actions 
of the bounce and of the original false vacuum; i.e., 



$ ^bounce Sj 

where the false vacuum horizon radius 



71 I A 



G 



h) > 



(3.12) 



(3.13) 



fv 



For actually seeking bounce solutions, it is convenient to recast our field equations a bit. 
First, the equations are simplified by defining 



Equation ( 13. 8j) then becomes 







1 d 



8nGdr 

If we define a mass A4(r) by 



K/-l)] + r 2 



this can be rewritten as 



f(r) = 
M' = 4vrr 2 



2GM(r) 



Ll\2 



Yf + v(<f>) 



(3.14) 

(3.15) 

(3.16) 
(3.17) 



2 These equations can also be obtained by substituting our ansatz for the field and metric into the full field 
equations. 
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Second, we note that substitution of Eq. (13.91) into Eq. (I3.10p gives 



= / 



r 



+ /y - 



dV 



(3.18) 



This and Eq. (I3.17|) give a pair of equations involving only / and (j). Once these have been 
solved, AB is readily obtained by integrating Eq. (13. 9p . With our convention that AB be 
equal to unity at the horizon, this gives 



A(r)B(r) = exp 



-8ttG / dr'r 



f 



(3.19) 



We must now determine the boundary conditions. We clearly need to avoid singularities 
at the center of the bounce, r = 0. This leads to the requirements 



/(0) 
f(0) 



1 





(3.20) 



The value of 0(0) is not predetermined, although it must be on the side of the barrier 
corresponding to the new vacuum (i.e., the true vacuum side for the case of false vacuum 
decay by bubble nucleation). 

The horizon is determined by the requirement that f(rn) = 0. Our prescription for 
specifying the temperature guarantees that the Euclidean metric is nonsingular at r#, even 
when the bounce is viewed as a topological four-sphere and not simply as the product of a 
three-ball and a circle. In order for <fi to also be nonsingular on the four-sphere we would 
have to require that 0'( r #) = 0. Equation (13 . 1 8[) shows that this implies that dV/d(f> = at 
the horizon, which in turn means that the field must be precisely at a vacuum value (or else 
exactly at the top of the barrier, as in the Hawking-Moss solution). For a nontrivial bounce 
this is only possible in the limit where the thin- wall approximation is exact. 

However, there is no reason to require such smoothness on the four-sphere. Instead, the 
only requirement is that following from Eq. (I3.18p . which implies that 







dV 



(3.21) 



r=ru 



if (j)' and <p" remain finite as r — > th- [In fact, the condition holds even for singular 0', 
provided that <fi' grows more slowly than (r# — r) 1//2 as the horizon is approached. Note that 
such a divergence does not imply a divergent action density, and in fact is consistent with a 
finite derivative of with respect to proper distance from the horizon.] 



IV. THE THIN- WALL APPROXIMATION 

The thin-wall approximation is often useful for gaining an intuitive understanding of 
tunneling problems. In this approximation, the bounce is approximated as a region of pure 
true vacuum surrounded by a region of pure false vacuum, with the variation of the scalar 
field confined to a thin wall separating the two. There are two essential requirements for this 
approximation to be valid. First, the thickness of the wall must be much less than the radius 
of the true vacuum region, so that the wall can be locally approximated as being planar. 
The second concerns the position dependence of the surface tension a of the wall. Although 
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a was independent of position in most previous tunneling calculations, position dependence 
is possible. However, for the thin-wall approximation to be reliable, the variation of a must 
be small on distance scales comparable to the width of the wall. 

In flat spacetime the thickness of the wall is determined primarily by the form of the 
potential barrier separating the true and false vacua. Reducing e, the difference in the 
energy densities of the two vacua, while keeping the shape of the potential barrier otherwise 
fixed increases the radius of the true vacuum region. By making e sufficiently small one can 
always reach the regime where the thin-wall approximation is valid. This is the case both for 
0(4)-symmetric zero-temperature bounces and for high-temperature solutions that describe 
a critical bubble in three dimensions. 

With gravity included, the situation is more complicated. Because the false vacuum hori- 
zon length limits the size of the bounce, the true vacuum region cannot be made arbitrarily 
large. To even have a possibility of the thin-wall approximation being valid, the shape of 
the potential barrier must be such that the natural width of the wall is small compared to 
the horizon length. If this is the case, then an appropriate reduction of e can ensure the 
validity of the approximation for the 0(4)-symmetric bounces. 

This turns out to not be true for our 0(3) xO(2)-symmetric bounces. To understand 
this, let us begin by assuming that the variations in the metric functions / and B as one 
passes through the wall are small enough that they can be treated as constants there. 
Equation (13.18P then reduces to 



If / is treated as a constant near the wall, this is identical with the analogous flat spacetime 
equation, except for a rescaling r — >■ f = r/yff. In other words, the wall thickness, measured 
in terms of r, is reduced by a factor of \ff compared to the flat spacetime case; this is to be 
expected, since it is f that measures the proper distance across the wall. Now let us further 
assume that the second, 2<f>'/r, term can be neglected. Multiplying the above equation by 
0' then leads to 



so that the quantity in brackets is constant through the wall. Evaluating this constant just 
outside the wall tells us that in the wall 




(4.1) 




(4.2) 




(4.3) 



Furthermore, we have 



wall 



A J w 



wall 



dr [V{(f>) - Vfy] = / 



ft** d(j) f 



r\2 




(4.4) 



where the integrals over r are restricted to the wall region. 
We can therefore define a surface tension 




(4.5) 
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and write the contribution to the action from the matter in the wall as 



Air /3 I drr 2 VAB 

J wall 



Vt 



fv 



(4.6) 



After the gravitational contribution to the wall action is included, Eq. (14. 6 j) agrees with the 
corresponding expression in Ref. j5|, verifying that the surface tension that we have defined 
in Eq. (14. 5 j) corresponds to that used in the latter paper. 

We must now go back and examine the validity of the assumptions that we have made. 
In particular, we must show that the fractional changes in r, /, and AB as one goes through 
the wall are all small. For the first of these, we note that the thickness of the wall is 



Ar 



<Atv ^2^(0)-^] 



(4.7) 



Just as in flat spacetime, this can be made arbitrarily small, without materially changing 
a, by making the potential barrier higher by a factor r y 2 > 1 and narrower by a factor I/7. 
Also, for fixed barrier height and width this expression decreases as one moves toward the 
horizon. 

Equation ( 13.191) shows that the fractional change in AB is 



A(AB) 
AB 



AirGr / dr f{<j>') 

J wall 



'\2 



2ttGW/. 



(4.8) 



If the mass scales in the potential are small compared to the Planck mass, this is also small. 
The difficulty is with the variation in /. From Eqs. ( 13 . 1 6[) and (I3.17P we find that 



A/ 



2G 



AM 



SirGr / dr 

J wall 



/(0f + v 



-8nGr aJf + V {v Ar 



Requiring this to be small compared to / gives the strongest constraint, 

8nGra < Jf < 1 . 



(4.9) 



(4.10) 



This will certainly fail if the bubble wall is too close to the horizon, where / vanishes. 3 
In particular, the limiting surface tension gn that is discussed in Ref. jBJ does not satisfy 
Eq. (I4.10p . In the field theory context, if the potential is such that Eq. (14. 5 p implies a }t a^, 
then the thin-wall approximation is not valid. 



V. LIMITING CASES 

Before we turn to our numerical results, it may be helpful to consider some limiting cases. 



3 On the other hand, if r is much smaller than m, Eq. (|2.8|) will be violated, in which case we expect to 
have multiple negative modes and therefore an unacceptable bounce. 
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A. Weak gravity 

Let us begin by focusing on the effect of varying G, with all other parameters held fixed. 
We first consider the case where G tends toward zero, so that gravitational effects would be 
expected to be small. If G is sufficiently small, <fi varies from a value near the true vacuum 
to one exponentially close to 0f v in an interval < r < f, where f is much less than Af v . 
In this region Eq. ( I3.18P is well approximated by the corresponding flat spacetime equation. 
Outside this region, integration of Eq. ( I3.17P gives 

M(r) = E m + ^^-, (5.1) 

where E& at is the energy of the critical bubble in flat spacetime. To leading order the horizon, 
where / vanishes, is at r# = Af v — GE^. The tunneling exponent, Eq. (13. 12ft . is given by 

B = 27rA fv £ flat = £ flat /T dS , (5-2) 

where T^s = l/27rAf v is the de Sitter temperature of the false vacuum. As expected, this is 
the Boltzmann exponent one would obtain for nucleation of a critical bubble in flat spacetime 
at a temperature T = Tds- However, because T <C 1/f, the high temperature bounce has 
multiple negative modes, and must therefore be discarded. 



B. Strong gravity 

Now consider the opposite limit, with G becoming large while all other parameters are 
held fixed. In this limit Af v decreases, leaving less and less space for the wall separating the 
two vacua. As Af v becomes smaller than the natural width of the wall, is restricted to an 
increasingly narrow range of values around the top of the barrier. Eventually, the solution 
becomes spatially homogeneous, with <p = top everywhere. This is, in fact, the usual 
Hawking-Moss solution, whose 0(5) symmetry has both 0(4) and 0(3)xO(2) symmetries 
as subgroups. 

The approach to the Hawking-Moss solution can be studied by taking 5(f> = (f) top — <p to be 
small. This quantity enters quadratically in the equations for the metric functions, Eqs. (13. 8 j) 
and (13.91) . so to leading order we can linearize the scalar field Eq. (I3.10P in 8<p and take the 
metric functions to be those of a de Sitter spacetime with V = V((p t o P )- The analysis is then 
quite analogous to the four- dimensional case that was considered in Ref . [l7[ , with the main 
difference being the replacement of Gegenbauer functions by Legendre functions. We find 
that the nontrivial bounce solution goes over to the Hawking-Moss solution when 4 

^^)<10. (5.3) 

-"top 



4 In [17] it was found that there is on 0(4)-symmetric solution merging with the Hawking-Moss solution 
whenever V" /H 2 at the top of the barrier is equal to N(N + 3) for any integer TV. For the 0(3) x 0(2)- 
symmetric solutions the corresponding critical values are 2 (TV + l)(2iV +1) — 2. 
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C. Large vacuum energy 

Another possibility is to keep G fixed, but to raise the vacuum energies by adding a 
large constant to V((f>). On a pure vacuum solution increasing the vacuum energy has the 
same effect on the metric as increasing G. We would expect the same to be true if the 
variation of V(4>) in the region between the two vacua is small compared to the absolute 
magnitude of the potential. As in the previous case, the bounce solution should go over to 
the Hawking-Moss solution when Eq. (15. 3p is satisfied. 

VI. NUMERICAL RESULTS 

To explore the properties of the bounces more closely, we considered a theory with a 
quartic potential with two minima. By shifting the zero of <p so that the minima are equally 
spaced about the origin, with tv = v and 0f v = —v, any such potential can be brought into 
the form 

V{(j>) = A (c v 4 - kv 3 <j) - ^v 2 <f) 2 + + ^ . (6.1) 

with < k < 1. The top of the barrier separating the two vacua is at = — kv. The theory 
is thus characterized by four dimensionless quantities: A, Co, k, and 

h = 8ttGv 2 . (6.2) 

The dependence on one of these, A, is rather simple. If 4>{x) = g(x) is a bounce solution for 
a given value of A, then <j)(x) = g(p(x) is a solution for the theory with A replaced by 7 2 A. The 
action of the bounce in the latter theory is 7 -2 times that of the original bounce. Thus, the 
classical solutions have a nontrivial dependence on only the three remaining dimensionless 
parameters. At the same time, it must be remembered that the validity of the semiclassical 
method requires that we be in the weak coupling regime A<1. 

We will find it useful to work with an alternative set of parameters, 

e = \{V fv -V ty ) = ^Xk 1 

V 4 6 

« = ^(V top -V iv ) = ±(3-k)(l + k) 3 , 

1 f 1 2k \ 

U = ^v = Af-- + y + C J , (6.3) 

that characterize the splitting of the vacuum energies, the height of the potential barrier, 
and the overall magnitude of the vacuum energy. In order to understand the dependence 
on these features of the potential, we will hold two of these (and h) fixed while varying the 
third. Such a variation corresponds to a more complicated path in the X-C^-k space, and 
may extend into the large A strong coupling regime. This does not matter for our purposes, 
since any such strong-coupling solution is related to a weak-coupling solution whose field 
and metric profiles differ from the original only by a rescaling of distance 

We have numerically solved the bounce equations for a variety of values of the parameters. 
Figure |5] shows the effect of varying G with the other parameters held fixed. The horizon 
radius decreases as gravity is made stronger, just as one would expect. For sufficiently 
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(c) 

FIG. 5: Evolution of the bounce as the gravitational constant is varied. The first two panels show 
as a function of (a) r and (b) r/m- In both cases the black dot indicates the value at the top 
of the barrier. The third panel shows / = 1/A as a function of r. Reading from left to right, 
the short-dashed purple, long-dashed brown, dotted red, dot-dashed blue, and solid black lines 
correspond to 8irGv 2 equal to 0.711, 0.704, 0.628, 0.251, and 0.126. In all cases U = 2, a = 3, and 
e = 0.3. 

strong gravity, the starting and ending points of the field move significantly away from their 
respective vacuum values. The solution eventually approaches the Hawking-Moss solution, 
with being constant and / precisely of the de Sitter form, in agreement with the arguments 
that we presented in Sec. [V] The Hawking-Moss solution is reached at the value predicted 
by Eq. 

We also argued in Sec.|V]that the effect of uniformly increasing Vf v should be qualitatively 
similar to that of increasing the strength of gravity. This is verified by our numerical 
solutions, which closely resemble those shown in shown in Fig. [51 Again, the Hawking-Moss 
solution is reached at the value predicted by Eq. (15. 3p . 

Figure [6] shows the effect of varying e, the dimensionless difference of the two vacuum 
energy densities. We see that decreasing e increases the bubble radius. This behavior is 
already familiar from the flat spacetime case, where the increased radius is needed for the 
outward pressure of the interior true vacuum to be able to overcome the surface tension of 
the bubble wall. Note, though, that we have included a curve for e = 0, corresponding to 
degenerate vacua, and that the outermost curve is for negative e; i.e., the case of a false 
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(c) 

FIG. 6: Behavior of the bounce as e is varied. Again, panels (a) and (b) show ^ as a function of r 
and of r/r/f , with the black dot indicating the value at the top of the barrier. Panel (c) shows the 
metric function / = 1/A. Reading from left to right, the dashed brown, dotted red, dot-dashed 
blue, and solid black lines correspond to e equal to 2, 1, 0, and —1. For all cases Uq = 2, a = 3, 
and 8irGv 2 = 0.126. 

vacuum bubble forming within a true vacuum region. Neither of these would be possible in 
flat spacetime, where the bubble radius tends to infinity as e is taken to zero, and where a 
true vacuum region cannot nucleate false vacuum bubbles. This is not so when gravitational 
effects are included. If the true and false vacua are both de Sitter, it is possible to "tunnel 
upward" from true vacuum to false [8|. Although the surface tension and the vacuum 
pressure would both tend to collapse the resulting bubble, the Hubble flow of de Sitter 
spacetime overcomes these if the initial bubble size is sufficiently large. 

In the CDL ansatz, the "new vacuum" region of the bounce is centered about the "North 
Pole" of the four-sphere, and the "old vacuum" region about the "South Pole" . Nucleation of 
a bubble of false vacuum rather than true is described by the same bounce solution, but with 
the labeling of the two poles simply interchanged. Matters are different for the 0(3) x 0(2)- 
symmetric bubbles. The spacelike slices of these are three-balls, with the new vacuum near 
the center and the old vacuum at the outer edge. There is no longer a symmetry between 
the two regions, and the bounce profiles for the two processes differ, as can be seen in the 
figure. 

Perhaps the most interesting behavior is that which occurs when we vary a, the parameter 
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FIG. 7: Variation of the bounce with barrier height. The panels on the left correspond to smaller 
values of a, with the short-dashed purple, long-dashed brown, dotted red, dot-dashed blue, and 
solid black lines corresponding to a equal to 1, 3, 5, 12, and 15. The right-hand panels correspond 
to larger values of a, with the short-dashed purple, long-dashed brown, dotted red, dot-dashed 
blue, and solid black lines corresponding to a equal to 15, 20, 25, 40, and 100. For all of these 
U = 2, e = 0.3, and 8irGv 2 = 0.126. As in the previous figures, the black dot indicates the value 
of the field at the top of the barrier. 
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characterizing the barrier height in the potential. Increasing the barrier height increases the 
surface tension in the bubble wall. In flat spacetime this, in turn, requires that the bubble 
nucleate with a larger radius. A second effect of the increased barrier height is that the 
bubble wall gets thinner. 

Figure H shows the effects of varying a on the 0(3)xO(2)-symmetric bounces. We see 
two distinct types of behavior. For smaller values of a, illustrated in the panels on the left, 
the scalar field starts close to the true vacuum value at r = and then comes exponentially 
close to the false vacuum by the time that the horizon is reached. In contrast with flat 
spacetime, the radius at nucleation decreases as a increases. This is largely due to the fact 
that the increased tension in the bubble wall causes the metric coefficient / to turn over 
sooner, leading to a smaller horizon radius. However, when the field profile is viewed as a 
function of r/r#, the result is closer to what the flat spacetime experience would have led us 
to expect. Whether viewed as a function of r or of r/r#, the bubble wall becomes thinner 
and steeper with increasing a. 

Viewing these results for small a, one might be led to expect that further increases in 
a would lead to an even steeper profile that pushed ever closer to the horizon, perhaps 
becoming a step function, with an infinitely thin bubble wall, for sufficiently large a. This 
is not what happens. Instead, the value of the field at the horizon, ^(rn), falls visibly short 
of the false vacuum by an amount that increases with a. [For the highest value of a shown 
here, 4>{th) = —.85?;.] Furthermore, the wall does not become thinner. Indeed, when plotted 
as functions of r/rjj, the various large-a field profiles are almost indistinguishable. 

VII. DISCUSSION AND CONCLUDING REMARKS 

We have studied a class of bounce solutions with 0(3) x 0(2) symmetry. These are curved 
spacetime analogues of the flat spacetime bounces describing thermal nucleation of critical 
bubbles. Following the interpretation of Ref . 0] , the critical bubble configuration is obtained 
by taking a constant r slice through the bounce. This slice is topologically a three-ball, and 
gives the scalar field and metric on a horizon volume. No reference is made to quantities 
beyond the horizon, and the end result of the process is the creation of one bubble, not two. 

Using a general quartic potential as an example, we have explored the behavior of the 
bounces in various regions of parameter space. Either increasing Newton's constant or raising 
V(<f>) by a uniform constant, with the shape of the potential otherwise held fixed, drives the 
bounce toward, and eventually merges it with, the Hawking-Moss solution. Decreasing the 
energy difference between the new true vacuum and the background false vacuum tends 
to increase the bubble radius at nucleation, but the effect is hardly as dramatic as in flat 
spacetime, even in the case where e changes sign so that the bounce is actually describing 
the nucleation of a false vacuum bubble in a region of true vacuum. 

The effect of increasing the surface tension in the bubble wall by raising the potential 
barrier is particularly notable. Initially, this makes the bubble wall thinner and moves it 
closer to the horizon. This is fits one's expectations from the behavior of the thin-wall 
analysis of Garriga and Megevand j3, [E| • However, with further increases the two analyses 
diverge. In the thin- wall/membrane case the bubble wall reaches the horizon at a critical 
surface tension a N . By contrast, in our field theory case the field profile at the bubble 
wall asymptotes toward one with nonzero thickness, and in fact shows little variation with 
increased surface tension. 

To be relevant for tunneling, it is not sufficient that the bounce be a solution of the 
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Euclidean field equations. In addition, the spectrum of fluctuations about the bounce must 
include one negative mode, to give an imaginary part to the energy of the false vacuum, 
but no more, since additional negative modes would signal that the bounce corresponded 
to a path that was not a local minimum of the tunneling exponent. There are potentially 
three types of negative modes for the 0(3) x 0(2) symmetric solutions. Two of these are 
closely analogous to the ones associated with thermal bounces in flat spacetime. The first is 
a mode that is independent of the Euclidean time and corresponds to the radial expansion 
or contraction of the critical bubble profile. We expect that this negative mode will always 
be present. The second class of modes are obtained by modulating the amplitude of the 
previous mode with a sinusoidal variation in the imaginary time. By analogy with the 
discussion in Sec. Ill we expect these additional modes to be present if the ratio of the 
bubble radius to the inverse temperature is too low. In order to avoid these, the bubble 
radius should be comparable to r H . 

The third type of negative mode is closely tied to the four-sphere topology of the bounce. 
It is perhaps easiest to visualize in the limit where the bubble radius is much less than the 
horizon radius. In this limit the shaded region in Fig. H] becomes a narrow band stretched 
along the circumference of the bounce. Moving this true vacuum region to either side reduces 
its length in the r direction and therefore reduces its action. This mode may be avoided if 
the bubble radius is much larger and the bubble wall is no longer thin but instead has a tail 
extending out toward the horizon. This is the same region of parameter space as is required 
to avoid additional negative modes of the second type. 

Finally, even if the 0(3)xO(2)-symmetric bounce has only a single negative mode, it is of 
practical relevance only if its action is less than that of the CDL bounce. Indeed, one of the 
goals of our investigation was to test the conjecture that the bounce with minimum action 
has 0(4) symmetry. It was already shown in Ref. jBf that in the thin- wall approximation 
the 0(3)xO(2)-symmetric bounce always has a higher action than the corresponding CDL 
action. Our numerical methods allow us explore the much larger parameter space that is 
allowed when one goes beyond this approximation. The results are summarized in Fig. |HJ 
Each of the four panels corresponds to a one-dimensional path through this parameter space. 
We see that the tunneling exponent for the 0(3) x 0(2) is never less than that of the CDL 
bounce. The two approach and become equal at large values of G and of Uq. This can be 
understood by recalling that the 0(3) xO(2)-symmetric bounce reduces to the Hawking- Moss 
solution when V" / H 2 < 10 at the top of the barrier [see Eq. (joTB"]) ]. For 4 < V" / H 2 < 10 
there is always a CDL bounce whose action is less than that of the Hawking- Moss solution, 
but at V"/H 2 = 4 the CDL solution merges into the Hawking-Moss. 

These numerical tests are, of course, not conclusive. Although we have no reason to 
expect it, there could be unexplored regions of our parameter space where the action is 
less than the CDL action. More generally, there may be other choices of the potential for 
which the CDL bounce is not optimal. Conclusive analytical results on this issue would be 
welcome. 
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FIG. 8: Comparison of the tunneling exponents for the CDL bounce (red diamonds) and the 
0(3)xO(2)-symmetric bounce (black circles) as one parameter is varied with the other three held 
fixed. When held fixed, U = 2, a = 3, e = 0.3, and 8irGv 2 = 0.126. 
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